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SOME PROBLEMS ON THE THEORY OF CREEP* 

By Y. N. Rabotnov 

Section 1. The term creep of metals is applied to the phenomenon in 
which, at temperatures beyond a certain limit, the metal subjected to a 
load slowly undergoes deformation with time. For the case of steel, the 
creep phenomenon must be taken into account at temperatures above 400° C. 
Very slow deformations for a prolonged period are cumulative and lead 
either to inadmissible changes in the dimensions of a structural part or 
to its failure. 

It is important to note that failure due to creep occurs for very 
small strains considerably less than those in static rupture. 

In the design of steam power units, boilers and turbines, creep is 
a basic factor which deter m ines the choice of the admissible stresses. 

On account of the extreme urgency of the problem, the creep phenomenon 
has claimed the widest attention of metallurgists, physicists, and to a 
lesser extent, technicians. 

At the same time, however, the theory of creep constitutes part of 
the mechanics of dense media and the mechanical formulation of the pro- 
blem may be given as the following: 

A body is subjected to the action of a given system of forces, or 
initial displacements are prescribed on its surface. It is required to 
find the stress distribution in the body and the changes of its defor- 
mations with time. 

Such a statement of the problem immediately raises the following 
question: What tests should be set up in order that the mechanical 

characteristics of creep (certain constants or functions) may be deter- 
mined? Is it sufficient for this purpose to make use of the generally 
accepted methods of testing, or is it necessary to supplement them? 

For the solution of the problem of creep as thus formulated, a 
mechanical theory of creep is required. Such theory, at the present 
state of knowledge of the physics of the process, must necessarily bear 
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an external, over-all character based, not on the investigation of 
microprocesses , but on the results of mechanical, tests. The existing 
physical theories are as yet far from providing a quantitative descrip- 
tion of the process in all its complication. 

The usual method of creep tests is that of obtaining the strains 
for a constant load. For small deformations it may be assumed that to a 
constant load there corresponds a constant stress. 

The results of tests are generally represented in the form of creep 
curves (fig. l), the time being laid off on the axis of abscissas and the 
strain e on the axis of ordinates . 

The intercept € Q represents elastic deformation if the stress d 

does not exceed the elastic limit of the material.. Often in place of 
the total strain there is laid off the plastic strain 

p = e - a / E 

Many attempts have been made to give an analytical expression for 
the creep curves . The different equations proposed may be divided into 
two groups: 


1. p = S(o)T(t) (l*l) 

2. p = g(o)e(t) + S(c)t (0(°°) = 1) (1.2) 

In writing equation (l.2), the essential assumption is that the rate of 
creep tends towards a constant value with time, that is, the creep curve 
has an asymptote. The function g(a) is the intercept e ' on figure 1. 
For the function S(a), the following equations have been proposed: 


S(a) = Aa n 

(Bailey) 

S(a) = ve°^ 

(Ludwik) 

S (a) = v(e a /^ - 1) 

(Soderberg) 

S(a) = 2v sh — 
M- 

(Nadai) 

S(a) = bae a /^ 

(Oding) 
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» For the function T(t ) , the best equation is apparently the power equa- 

tion 


T(t) = t m (0 < m < 1) 

If equations of type (l.2) are used, the function 0(t) assumes 
the following forms: 

e(t) = 1 - e _pt (Mac Vetty) 

0 (t) = a ' V ' t (Oding) 

For g(cr), it has not as yet been possible to establish a law. 

The series of creep curves for different values of a give the 
representation of the functional relations between the three variables 
o, e , and t. If € and a are laid off on the coordinate axes, a 
series of curves is obtained, which is shown in figure 2 and character- 
ized by different values of t. Figure 2 refers to the test data of 
Robinson which are unique in that these tests extended for 100,000 hours 
(from March 27, 1931 to October 8, 1942). We worked over a very large 
number of the tests in the series by this method and obtained in 1 
cases, with an accuracy not exceeding the limits of the experimental 
accuracy, affinely related curves in the ae plane. 

On the basis of this result, the following formula describing the 
law of creep for constant stress is proposed: 

<¥>(0 = [l + a (1.3) 

For the function G(t), good results are generally given by the following 
expression: 

G(t) - t 1 - 11 

where the coefficient a fluctuates about the value 0.7, while the 
coefficient X changes very strongly for different materials . 

The scatter in transferring all the points on to the curve t = 0 
(fig. 2) is found to be not greater than the scatter of the modulus of 
elasticity of the specimens on which the given series of tests was con- 
ducted. 
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Section 2. The mechanical theory of creep for the case of a single 
axis must establish such relation between e and a, containing the 
time or time operators, which would permit predicting the course of the 
process varying with time. In particular, a second extreme case of the 
one-dimensional problem is the. problem of relaxation. The latter is the 
process of decrease in stress in a rod the length of which remains con- 
stant. Test data on relaxation sire very meager, only those data being 
of value which, together with the curves of relaxation, give the curves 
of creep for the same material. The only reliable data in the general 
literature axe those published by Davis in 1943 on copper. 

The various mechanical theories of creep existing at the present 
time may be divided into three groups: 

(l) Theory of constant rate: Assuming the existence of an asymp- 

totic curve of creep, the curve is replaced by a straight line parallel 
to the asymptote and intercepting the segment € Q = c/e on the axis of 

ordinates . Then 


P = S(c) (2.1) 

The theory of constant rate assumes this as the true relation for any 
conditions. In particular, for € = constant, p = -&/E, and from equa- 
tion (2.1) there is readily obtained the law of relaxation which grossly 
contradicts test data, since the neglect of the primary creep (the curvi- 
linear part of the curve £t) is not permissible in the problem of 
relaxation. In the case of the nonuniform stress state, equation (2.1) 
leads to very great difficulties and even the problem of the pure bending 
of a rod of rectangular cross section is not solvable. 

The greater number of authors employing this theory take still 
another step and neglect the elastic deformation. The fundamental equa- 
tion is then the following: 


£ = S ( 0 ) (2.2) 

Equation (2.2) is generally put at the basis of the theory of sec- 
ondary creep widely applied in technical computations . Below shall be 
given another basis for this theory which considerably generalizes it. 

For the present we may note that from the point of view of equation (2.2) 
the problem of relaxation has no significance. 

(2) Theory of aging: The theory of aging postulates the existence 

of a definite relation between a, €, and t: 
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f(a,€, t) = 0 (2.3) 

or between p, a, t 

F(p, a, t) = 0 (2.4) 

It is easily 6hown that any theory, the fundamental equation of which 
contains the time explicitly, is contradictory. The physical law must 
be invariant relative to a time origin. Applying the theory of aging 
to the successive loading and unloading may yield absurd results. 

The theory of aging in the form of equation (2.3), however, gives 
for smoothly varying loads satisfactory agreement with experiment. The 
methods of computation based on it are relatively simple and at the same 
time permit taking into account all the characteristic experimental 
curves which may be obtained in tests. Hence, one of the variants of 
the aging theory can be recommended as a technical method of computing 
structural parts working under the conditions of high temperatures. This 
point of view has been developed by us in a paper presented at the 
session of the Soviet Academy of Sciences in March 1948. 

(3) Theory of strain hardening: This theory postulates the exist- 

ence of an unvarying relation among the rate of plastic deformation, its 
magnitude, and the magnitude of the stress: 


*(p> P> o) = 0 (2.5) 

Methods exist for the graphical construction of the relaxation curve 
by a given family of creep curves on the basis of hypothesis (2.5). The 
analytical formulation of this theory, the choice of the functional 
relations for which it is possible to integrate equation (2.5), is found 
to be very difficult. 

Section 3. The theories of creep enumerated in section 2 axe not 
capable of explaining a number of phenomena observed during experiment. 
There is first of all the case for the strain- hardening effect. A 
specimen initially strained by a large force evidences creep to a con- 
siderably less extent than a specimen not initially strained. If the 
specimen was tested under stress cj^ and is decreased to 02 , the 

creep practically vanishes. If, however, the stress 0% was initially 

imposed, the creep for this stress may be very marked. There are no 
serious investigations of this problem in the literature. The published 
data undoubtedly give a qualitative account of the phenomenon but the 
quantitative aspect still awaits investigation. 
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The second effect is the so-called reverse creep. If a specimen is . 

subjected to a creep test at constant load and the load is then removed, 
the specimen immediately shortens by the amount of the elastic elongation 
but the process does not stop there. In the course of time the specimen » 

continues to shorten, returning in this way a part of its residual 
deformation. This phenomenon of the type of elastic aftereffect has been 
subjected to a careful experimental study but not one of the previously 
enumerated theories provides an explanation for it. 

The theory of creep proposed by us represents an extension of the 
theory of elastic heredity of Volterra to plastic deformation. In the 
same way as the elastic heredity develops about a straight line a = Ee 
in the cr € plane, the plastic heredity, or creep, develops about a 
certain curve in this plane. It will be inconvenient, in what follows, 
to lay off as usual the value of the stress a on the ordinate axis and 
we shall therefore take a certain fictituous plane cpe and a curve 
cp = cp(e) in this plane (fig. 4). 

The curve op(€ ) represents the ideal curve of strain with the 
exclusion of the time factor (actually never realized). For active 
processes, that is, those accompanied by motion along this curve upward, 
the fundamental law is written in the following manner: 


cp(c) = (1 + K*) a (3.1) 

where K* is the integral operator of Volterra, that is, 

K*a = K(t - T) o(t) dT 

J o 

Here and in what follows, use will be made of the notations and 
results of our previous paper (ref. l). 

The increase in the stress a is no longer a criterion of the 
activity of the process as in the theory of plasticity but such criterion 
is given by the increase in e. 

For unloading processes, it is necessary in the left-hand side of 
equation (3.1) to introduce in place of cp(c ) the value of the ordinate 
of the linear unloading AB expressed as a function of € : 


E (c - €•) + cp' = (1 + K *) a 


(3.2) 
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If r* is the solving operator, 

j: _ i _ r * 

1 + K* 

there follows from equation (3.1) 


a = (1 - r*) cp(c) 


Setting 


K* • 1 = G(t) and T* * 1 = R(t) 
the law of creep will then he 


cp(e) =[l + G(t) 3 ° 0 (3.3) 

and the law of relaxation 


a = [1 - R(t)] cp(e 0 ) (3.4) 

Equation (3.3) entirely agrees with that obtained from tests, for- 
mula (l.3). The expression given in section 1 for G(t) shows that the 
operator K* possesses a singularity of the type of the Abel operator. 
In the simplest case we may assume 


k* =■ 


where I*_ a is an operator with kernel (t - T )" a /T(l-a) . Introducing 

the 3 -operators employed in the previously cited paper (ref. l) yields 
a very close agreement with test results, since the 3 -operator contains 
an additional constant. [HACA Reviewer's Rote: The 3 -operator is an 

“ ft n/, %n(l+a) 

operator with the kernel 3 a (p, t - t) = (t - t) £_/ f [ ( n ~+ 1 ) (1 + a )] * 
The additional constant referred to is |3.] 

If, as is usually the actual case, the body is acted upon by con- 
stant loads or the displacement of its points is maintained by constant 
stresses , it is necessary in actual computation to deal with the two 
functions G(t) and R(t) which may be determined from tests forgetting 
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about their origin from the kernel and as resolvents of the same inte- 
gral equation. 

This variant of the technical theory of creep is a forward step as 
compared with the aging theory which retains equation (3.3) both for 
creep and for relaxation. 

The phenomenon of reverse creep (fig. 5) will now be considered. A 
stress o is applied to the specimen and it undergoes the instantaneous 
deformation o/E (point A)j in the course of time T, creep ocpurs along 
the curve ABj at the instant T the stress is removed and at the same 
time there is removed the elastic deformation BC = OA = a/ E. In the 
time elapsed 0 after the unloading there is also removed the defor- 
mation € r (point D) . This process will be followed in the plane cpe 

(fig. 6). 

The instantaneously applied load corresponds to the motion along 
the curve from point 0 to point A; during time T the strain increases 
and the function cp(e) up to point Bj the creep process is described 
by equation (3.1). At point B 


cp(e B ) * [l + G(T)] a (3.5) 

With instantaneous unloading we drop to point C on the segment 
BC = AO, the reverse creep corresponding to the motion along this line 
up to the point D at instant 0. [ MCA Reviewer's Rote: It is also of 

interest to examine the effect of an instantaneous increase in load 
occurring after a period of constant- stress creep. If e Q is the strain 

at the end of the creep period and 5e is an instantaneously imposed 
strain increment, then, according to equation (3.1), there is an 
instantaneous stress increment given by 


bo = cp( € 0 + 8e ) - cp(e Q ) (a) 

This result may be at variance with the facts, for recent experimental 
evidence, including tests on the propagation of plastic waves in bars 
subject to creep, suggests that after some creep has taken place materials 
behave elastically for small instantaneous increments of stress, that is, 

00 
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Only for sufficiently small € 0 would equation (a) reduce to equa- 
tion (b).] Here equation (3.2) must be used: 

E (c - e B ) + cpg = [G(T + 0) - G(0)] a (3*6) 

[MCA Reviewer's Note: The coefficient o in equation (3.6) was 

incorrectly omitted in the original.] 

It is seen from figure 6 that 


c r = 



Making use of equations 


(3.5) and (3.6) yields 


«r 



+ G(e) 



(3.7) 


The linear dependence ( of the reverse creep on the stress and the 
symmetry of its dependence on T and 0 is well confirmed by experi- 
ment (ref. 2). 


An important particular case 
cp(e) = a«E and K* = XI* „. The 

QL —CL 


of the relation (3.1) 
law of creep is then 


is obtained when 


ae^ = 
a 


It. X 4-1 ” C. 

L 1 + * 


For large t, the second term in the brackets is the dominating one and 


e = Acr^t' 


m 


(3.8) 


where 


1 m 1-a 

n = P ' m = ~ 


rate. 


If 1 - a = p, formula (3.8) gives, in the limit, a constant creep 


Section 4. The application of the theory presented in the preceding 
section to the problem of pure bending is considered herein. For sim- 
plicity, we restrict ourselves to the case of a rod of rectangular cross 
section, since we are concerned with the theoretical side of the problem. 
Assuming the hypothesis of plane sections, set 
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e = z/p 

where p is the radius of curvature of the bent axis, and z the 
coordinate measured from the neutral axis and varying from -h to +h. 

By equation (3.l) 


(1 + K*) a = cpf-\ 

\ p / 


(4.1) 


Multiplying by bzdz, where b is the width of the section, we 
integrate from z = 0 to z = h, and multiply the result by 2 (the 
function cp(e ) is analytically determined only for positive values of 
the argument; it is continued as an odd function in the region of nega- 
tive values). Since the time and space operators are interchangeable, 

^h 

(1 + K*) M = 2b cpf-j zdz 
Jo ' ' 

where M is the bending moment. Introducing the notation 


1 

x2 


rx 


j 


o 


cp(x) xdx = cp(x) 


yields 


(1 + K*) M = 2bh 2 cp(tj (4.2) 

The graph for the function cp(h/p) can be easily constructed. It 
is thus always possible to find the magnitude of the curvature from the 
given moment M for a given instant of time and then, by solving inte- 
gral equation (4.1), to determine the stress distribution. 

If the function cp(e) is a power function, the problem is consi- 
derably simplified. Let cp(x) = ax^ . Then 


cp(x) 


2 + P 


Dividing equation (4.l) by (4.2) yields 
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1 + K*) o 2 + (3 /z\ P 

1 + K*) M 2bh 2 \h) 


whence 


a 


2 + p 
2bh 2 



M 


(4.3) 


The distribution of the stresses is found to be independent of the time. 
We thus have in a certain sense a secondary creep although the rate of 
deformation is not constant as is clear from equation (4.2). 

The result, equation (4.3), was obtained by Davis (ref. 3) from 
entirely different considerations. He made use of the theory of strain 
hardening, using equation (2.4) in the form 


• TV_k t 

p = rip cf 

As he found difficulty in accurately solving the problem as thus 
posed, Davis identified the plastic deformation with the total defor- 
mation and arrived at a creep law of the type (3.8) from which he readily 
obtained equation (4.3) in the case of bending. The experimental verifi- 
cation conducted by Davis well confirms the theoretical result. 

Section 5. The generalization of any of the theories of creep to 
the three-dimensional case may be effected if there are taken as a basis 
the equations of the theory of small elasto- plastic deformations: 

2a. a. 

ax “ a 3€~ Cx ' T xy = 3e b € *y ( 5 * 1 ) 

The change in volume at high temperatures may be neglected with 
greater justification than for normal temperatures since, starting from 
400°, the Poisson coefficient for steel assumes a value of the order of 
0.45 to 0.47. 

The intensity of the stresses and the strains which are determined, 
following A. A. Ilyushin, as 
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£ i 



o z ) 2 + (o z - ° x ) 2 + 6 (r 2 y + T yz + 4) 


C z) 2 + ( £ Z - £ X ) 2 + | ( £ ly + 4 + 4) 


are assumed to be connected with each other by the same relations as the 
stress and strain in the one-dimensional problem. 

By the theory herein. 


cp(Ci) - (1 + K*) Oi (5.2) 

The problem of a pipe under the action of an internal pressure is 
very simply solved on the assumption of the absence of axial deformation. 
Let a r and Oq be the radial and transverse stresses, £ r = du/dr 

and €q = u/r the corresponding strains expressed in terms of the radial 

displacement u, b and a the external and the internal radius of the 
pipe, and q the internal applied pressure. Since C x = 0, there 

follows from the condition of incompressibility 


whence 


du u 
dr + r 


0 


& 





where e is an as yet undetermined function of time. Equations (5.1) 
give [NACA Reviewer's Note: In the following three equations, an error 

in sign appearing in the original has been corrected.] 
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whence 


Cfp ” O r\ — 


“VS 01 

Substituting (5.3) in the equation of equilibrium 


(5.3) 


da a - a 

+ -Z- = 0 

dr r 


yields 


da _ a. 

r 2 i 

~V3 r 

Multiplying by (l + K*), integrating from r = a, and noting that 


(5.4) 


yield 


Setting 


(1 + K*) a ± -qp(Ci) 


n.r 


(1 + K*)(a r + q) - 


,/3 


a 


ea^\ dr 


W ~ 


Qp(x) ^ = $>(x) 


gives 


(1 + K*)(a r + q ) =jz 

V 3 


9(e) - 


(5.5) 


For r = b, a r = 0. From this condition, the integral equation for 
determining the function of time e is obtained: 


(1 + K*)q = jz 

Vs 


$(e) - $^e 


( 5 . 6 ) 
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The case where cp(x) = axP here too leads to considerable simplification; 
in this case 


% 

$(x) = | xP 


[NACA Reviewer's Note: The parameter a in axP is evidently being 

taken as 1.3 

From equation (5.5), 


(1 + K*)(a r + q) 



a 2 * 

r 2(B 


(1 + K*)q 



a 2 P~ 

b 2 P 


Dividing one equation by the other cancels the time function e on 
the right side. The integral operator (l + K*) likewise cancels and a 
time- independent law for the stress distribution is obtained: 


/a\ 2|J r 22 - b 22 
- l(r) b 2|3 . a 20 


(5.7) 


The velocities are not constant and they are found from equation (5.6) 
where, if q is given, the case reduces to quadratures and the succes- 
sive determination of the function e from the graph. If e is given, 
the problem reduces to the solution of an integral equation for which it 
is necessary to know the resolvent of the kernel. 


Section 6. Proceeding to the general case of the three-dimensional 
problem, we restrict ourselves to the consideration of those processes 
in which the stresses vary simultaneously in proportion to a parameter 
X(t), while the strains vary in proportion to a parameter p(t). Then 


°X = A °x 


€ v = 


where a x , . . . e" . . . are functions only of the coordinates. 

[NACA Reviewer's Note: The bars which appear above o x and e x in th: 
clause were incorrectly omitted in the original.] Equations (5.1) give 

7 

v 

M°x 


— N 2 x 

' a ) = 3 a i T. 

1 


£S~Z 
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The stresses ck axe connected with the strains by the relation (5.2). 
If cp(€.,) is a power function 

= a 4 = an p F| (6.2) 

Multiplying both sides of equation (6.1) by 1 + K* gives 

(1 + K*) A (cr x - d) = \ |/ a 7^ ~ (6«3) 

3 € i 

[NACA Reviewer's Note: Errors in the subscripts on the right side of 

equation (6.3) appearing in the original have been corrected.] 

Equations (6.3), of which only the first is written out, correspond to 
the integral equation 


(1 + K*) A = (6.4) 

and to a system of relations entirely agreeing with the equations of the 
theory of small elasto- plastic deformation: 





(6.5) 


where 


°i = a ei = cp(€ i ) (6.6) 

Systems (6.5) and (6.6) are systems of equations of the theory of 
secondary creep taken in a more general sense than the ordinary. The 
rates are now no longer constant but variable since the factors A and 
M- which depend on the time are obtained from the integral equation (6.4). 
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